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Abstract 

Fermiophobic Higgs bosons (hf) exhibiting large branching ratios to two photons 
can arise in models with two or more scalar doublets and/or triplets. In such 
models the conventional production mechanisms at hadron colliders, which rely 
on the hfVV coupling {V = W,Z), may be rendered ineffective due to severe 
mixing angle suppression. In this scenario, double hf production may proceed via 
the complementary mechanism qq' —> H±hf with subsequent decay —> hfW*, 

leading to events with up to 4 photons. We perform a simulation of the detection 
prospects of hf in the multi-photon (> 3) channel at the Fermilab Tevatron and 
show that a sizeable region of the (m H ±, rrih f ) parameter space can be probed during 
Run II. 
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1 Introduction 


Neutral Higgs bosons with very suppressed couplings to fermions - “fermiophobic Higgs 
bosons” (hf) [I]- may arise in specific versions of the Two Higgs Doublet Model (2HDM) 
or in models with Higgs triplets |I]. Such a hf would decay dominantly to two 
photons, hf —> 77 , for m/ l/ < 95 GeV or to two massive gauge bosons, hf —> VV^*\ 
( V = W ± , Z ) for rn hf > 95 GeV 0 El- The large branching ratio (BR) for h/ —> 77 
would provide a very clear experimental signature, and observation of such a particle would 
strongly constrain the possible choices of the underlying Higgs sector u 00 la 13 nni nu- 

Experimental searches for hf at LEP and the Fermilab Tevatron have been negative 
so far. Mass limits have been set in a benchmark model which assumes that the coupling 
hfVV has the same strength as the Standard Model (SM) Higgs coupling VV<f)°, and that 
all fermion BRs are exactly zero. Lower bounds of the order rrih f > 100 GeV have been 
obtained by the LEP collaborations OPAL[T2|. DELPHI|13|. ALEPH^;, and L3jl51. 
utilizing the channel e + e _ —> hfZ, hf —> 77. At the Tevatron Run I, the limits on 
rrih f from the D0 and CDF collaborations are respectively 78.5 GeV |16| and 82 GeV 
na at 95% C.L., using the mechanism qq' —> V* —> hfV, hf —> 77, with the dominant 
contribution coming from V = W ± . For an integrated luminosity of 2 fb 1 , Run II 
will extend the coverage of rrif- LJ in the benchmark model slightly beyond that of LEP 
IEH2DI. hr addition, Run II will be sensitive to the region 110 GeV < m/ l/ < 160 GeV 
and BR (hf —» 77) > 4% which could not be probed at LEP. A preliminary search in the 
inclusive 2y + X channel has been performed with 190 pb^ 1 of Run II data eh 

However, the hfVV coupling in a specific model could be suppressed relative to the 
(j)°VV coupling by a mixing angle, leading to a weakening of the above mass limits. If this 
suppression were quite severe (hfVV/(j)°VV < 0.1) a very light hf ( rrih f << 100 GeV) 
would have eluded the searches at LEP and the Tevatron Run I in production mechanisms 
which rely upon the hfVV coupling. Therefore it is of interest to consider other production 
mechanisms for hf which may allow observable rates if the hfVV coupling is suppressed. 
Since the couplings hfVV and hfVH (where H is another Higgs boson in the model) are 
complementary, two LEP collaborations, i.e. OPAL (QJ and DELPHI [121, a % () searched 
for fermiophobic Higgs bosons in the channel e + e _ —> A°hf, and ruled out the region 
rriA + mh f < 160 GeV. However, a very light rrih f < 50 GeV is still possible if tyia is 
sufficiently heavy. 

An alternative production mechanism which also depends on the complementary 
hfVH coupling is the process qq' —> H ± hf [221, |2Sj- Such a mechanism is exclusive 
to a hadron collider, and can offer promising rates at the Tevatron Run II provided that 
H± is not too far above its present mass bound m#± > 90 GeV. This alternative ex¬ 
perimental signature depends on the decays of H ± . I 11 fermiophobic models the decay 
R ± -»■ hfWM can have a larger BR than the conventional decays H ± —> tb, tv mm, 
which leads to double hf production. 

In this paper we analyze the inclusive production of multi-photon ( 37 ’s or 4y’s) final 
states at the Tevatron RUN II via the mechanism: 

pp —> hf W =b —> hfhfW ± —* 777 ( 7 ) + A" . 
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In the 2HDM the multi-photon signature arises in the parameter space rrih f < 90 GeV, 
rriH± < 200 GeV, and tan/? > 1. In this region, BR(hf —■> 77 ) ~ 1 and BR( y H ± —> 
hfW * ± ) ~ 1, leading to a 47 + leptons or jets signature. The multi-photon signature has 
the added virtue of being extremely clean concerning the background contamination, in 
contrast to the conventional searches for single hf production in the channels 77 + V and 
77 + A". In the present work we show that the multi-photon signal can be observed in 
a large fraction of the rrif lf 8 ) m#± plane at the Tevatron RUN II. I 11 fact, at 3a level of 
statistical significance, the RUN II will be able to exclude Higgs masses up to m H ± <240 
GeV for very light m ^., or rrih f < 100 GeV for m H ± ~ 100 GeV. 

Our work is organized as follows. In Section 2 we give a brief introduction to fermio- 
phobic Higgs bosons, exhibiting the main decay channels of hf and H ± . The possible 
fermiophobic Higgs production mechanisms and respective signatures are described in 
Section 3. We present our analyses in Section 4 and Section 5 contains our conclusions. 


2 Fermiophobic Higgs bosons 

In this section we briefly review the properties of hf. For a detailed introduction we refer 
the reader to 0,0, mu, 'Hi- Fermiophobia can arise in i) 2HDM (Model I) and ii) Higgs 
triplets models. In (i) the imposition of a discrete symmetry together with a vanishing 
mixing angle ensures exact fermiophobia at tree-level. In (ii), gauge invariance forbids 
any coupling of hf to quarks while lepton couplings are strongly constrained by neu¬ 
trino oscillation data and lepton flavour violation experiments, resulting in approximate 
fermiophobia at tree-level. 


2.1 2HDM (Model I) 


If $1 and $2 are two Higgs SU( 2) doublets with hypercharge Y = 1, the most general 
SU( 2) x U( 1) gauge invariant scalar potential is |2By 

V = 1 + m 22 $ 2$2 - (mi 2 4>{T 2 + /i.c.) + |Ai 

( 1 ) 


+|a 2 (^$ 2 ) + a 3 ($ 1 $!) (^$ 2 ) + a 4 ($I$ 2 ) 

+ | TJA 5 ^ < f > | < f ) 2 ^ + Ag + A 7 (^ < I ) 2 < I ) 2 ^ d ) I < I ) 2 + h.c 


If the discrete symmetry $1 —> —$1 is imposed the couplings Ag = A 7 = 0. However, the 
term proportional to m\ 2 can remain as a soft violation of the above discrete symmetry 
and still ensure that Higgs-mediated tree-level flavour changing neutral currents are absent 
j3j . Note that the above 2HDM potential contains one more free parameter than those 
studied in Refs. mm- 

The potential in eq. 0 breaks SU{2)l x U(1)y down to U{l) em when the two Higgs 
doublets acquire vacuum expectation values 


<*i> 


1 

71 



($ 2 ) 


1 

71 



( 2 ) 
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which must satisfy the experimental constraint m 2 z = |( g 2 + g' 2 ){y\ + v\ ) ~ (91 GeV) 2 . 
The minimization conditions that define the vacuum expectation values in terms of the 
parameters of the potential (A 6 = A 7 = 0) are 


ti = m 2 u v i - ml 2 v 2 + |Aiuf + |(A 3 + A 4 + \ 5 )v 1 v% = 0 
t- 2 = m\ 2 v 2 - m\ 2 v i + \X 2 v\ + |(A 3 + A 4 + A 5 )u 2 w 2 = 0 


( 3 ) 


from which m\ x and m\ 2 can be solved in favour of m 2 z and tan/? = V 2 /V 1 . 

The neutral CP-odd Higgs mass matrix is, after using the minimization conditions, 


2 = ( m 2 12 tp - A 5 v 2 s 2 0 ~rn\ 2 + X 5 v 2 spcp A 

-m\ 2 + A ->v 2 spcp rri\ 2 /tf 3 - A 5 v 2 c 2 p ) 


m = 


( 4 ) 


and is diagonalized by a rotation in an angle /?. We define sp = sin/3, cp = cos/3, and 
tp = tan/3. M 2 4 has a zero eigenvalue corresponding to the neutral Goldstone boson 
while its second eigenvalue is the mass of the physical CP-odd Higgs boson A, 


m 


m A = 


12 


SpCp 


X 5 v 2 


( 5 ) 


with v 2 = v 2 + v 2 . The charged Higgs mass matrix is given by 


M 2 h± = 


_ ( rnj 2 tp - |(A 4 + A 5 )v 2 s 2 p -m\ 2 + |(A 4 + A 5 )u 2 s^\ 

-rn\ 2 + |(A 4 + A 5 )v 2 spcp m\ 2 /tp - |(A 4 + X 5 )v 2 c 2 p ) 


( 6 ) 


which also is diagonalized by a rotation in an angle /?. It has a zero eigenvalue corre¬ 
sponding to the charged Goldstone boson, and the charged Higgs mass is 


m 2 H ± = m 2 A + -( A 5 - A 4 )u 2 


( 7 ) 


Here we see that the charged and the CP-odd Higgs masses are independent parameters, 
as opposed to supersymmetry, where the mass squared difference is equal to m ^ at tree 
level. 

The neutral CP-even Higgs mass matrix is given by 


M^ 0 = 


—m\spcp + (A 3 + A 4 )u 2 s / 3 c /3 \ 
-m\spcp + (A 3 + X A )v 2 spc-p 7n 2 A c 2 p + X 2 v 2 s 2 p 1 


m A si + Ai vcp 

,2 n „ , /\ 1 \ \„,2, 




( 8 ) 


and the two eigenvalues are the masses of the neutral CP-even Higgs bosons h and H. It 
is diagonalized by an angle a defined by 


sin 2 a = 


[— m vl + (A 3 + A 4 )u 2 ] S2P 


m 2 A C 2 S - A, iPcp + \ 2 v 2 s 2 0 + [m 2 A - (A 3 + A 4 )w 2 ] 2 s 20 


( 9 ) 


Fermiophobia is caused by imposing the mentioned discrete symmetry $1 —> —$1 
which forbids $1 coupling to the fermions. This model is usually called “Type I” |2j. 
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However, fermiophobia is erased due to the mixing in the CP-even neutral Higgs mass 
matrix, which is diagonalized by the mixing angle a, when both CP-even eigenstates h° 
and H° acquire a coupling to the fermions. 

The fermionic couplings of the lightest CP-even Higgs h° take the form h°ff ~ 
cos a/ sin/3, where / is any fermion. Small values of cos a would strongly suppress the 
fermionic couplings, and in the limit cos a —> 0 the coupling h°ff would vanish at tree- 
level, giving rise to fermiophobia. This is achieved if 

m\ = (A 3 + A 4 )u 2 . (10) 

Despite this extra constraint, the parameters rriA, rnH±, and tan/3 are still independent 
parameters in this model. 

However, at the one-loop level, Higgs boson couplings to fermions receive contributions 
from loops involving vector bosons and other Higgs bosons, 

/ 

~ )™- f C 0 (ml, 0, 0; 0, m 2 w , m 2 w ) 

1 

where we have naively estimated the contribution of the loop with the Passarino-Veltman 
function Cq. To get an order of magnitude of the correction we approximate Cq ~ 1 /m\, 
expected in the limit of large Higgs mass, and compare this correction with the tree-level 
vertex in the SM g<pOff ~ gmj/2mw ■ We find 

HD 

g<ffi f f 647T 2 V m h ) 

This estimation is also applicable if mh ~ wiir replacing rrih by rriw- This is a very small 
correction. Nevertheless, we note that the proper renormalization of the 0°// vertex 
involves a counterterm that has to be taken into account. It is conventional to define an 
extreme hf in which all BRs to fermions are set to zero. This gives rise to benchmark 
BRs which are used in the current searches to set limits on mh r 

2.2 Higgs Triplet Models 

Fermiophobia (or partial fermiophobia) can arise for scalar fields in isospin 1 = 1 triplet 
representations. Gauge invariance forbids any couplings of the triplet fields (y) to quarks. 
For hypercharge Y = 2 triplets, the neutral Higgs held y° can couple to leptons (vi 7) via 
the following Yukawa type interaction M - 1 

+ h.c . (12) 

1 Note that there is no such interaction for Y = 0 triplets, which are rendered fermiophobic as a 

consequence of gauge invariance 
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Here h l3 {i, j = 1, 2, 3) is an arbitrary coupling, C is the Dirac charge conjugation operator, 
ipn = (zy, 4)1 is a left-handed lepton doublet, and A is a 2 x 2 representation of the Y — 2 
complex triplet helds: 


A 


X + /V2 X ++ \ 

X° -X + /V2 ) ' 


(13) 


The interaction described in eq. m has the virtue of being able to provide neutrino 
masses and mixings consistent with current neutrino oscillation data, without invoking a 
right-handed neutrino. If the real part of the neutral triplet held y 0r acquires a vacuum 
expectation value (vev) (y 0r ) = b , the following Majorana mass matrix (m y ) for neutrinos 
is generated: 

rriij = V2h i: jb . (14) 


Neutrino oscillation data constrain the product hijb , while hij is constrained directly by 
lepton flavour violating processes involving p and r e.g. fi —> ey, /j —> eee [2E,. Hence it 
is clear that y° is partially fermiophobic, with a small coupling to neutrinos. 

We will consider the Higgs Triplet Model (HTM) of reference [3| in which a complex 
Y = 2 triplet (A) and a real Y = 0 triplet (£ + ,£°,£ - ) are added to the SM Lagrangian. 
The HTM preserves p = 1 at tree-level if the vev’s of both the neutral members are equal 
(x°) = (£°) = b. Taking the vev of the Higgs doublet (<h°) = a, one has the following 
expression for mw 

m w = \d 2 (a 2 + 8b 2 ) = i g 2 v 2 (15) 


where v 2 = 246 2 GeV 2 . It is convenient to define a doublet-triplet mixing angle analogous 
to tan 0 in the 2HDM 


sindf/ = 


8 b 2 

a 2 + 8 b 2 


1/2 


(16) 


In the HTM the physical Higgs boson mass spectrum is as follows (in the notation of 

) 

rr±± rr± tt 0 rr± ttO rrO rrO' (T7\ 

-^5 i^5 T5>^3 )^3TlTl • GO 


The first five scalars are mass eigenstates, while the latter two can mix in general; see 
below. Hi plays the role of the SM Higgs boson and is composed of the real part of the 
neutral doublet held. The eigenstate is entirely composed of triplet helds and is given 
by 

ff?' = -k(^ x ° r +,e°). (is) 

From the theoretical point of view, the size of the triplet vev b is only constrained by 
the requirement that the doublet vev a is sufficiently large to allow a perturbative top 
quark Yukawa coupling. However, experimental constraints on sin 6h can be obtained by 
considering the effect of Hf on processes such as b —> sy, Z —> bb and B — B mixing jTT)! . 
Since Hf has identical fermionic couplings to that of H ± in the 2HDM (Model I) with 
the replacement cot /3 —>• tan Oh, one can derive the bound sin 0h < 0.4. 
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In eq. m the x° r component in Hf couples to vv via the coupling. One can 
see from eq. m that the decay H® — > 77 , mediated by W loops proportional to b, will 
dominate over H ° > vv if b is of the order of a few GeV. Thus H® is a candidate for a 
hf, with BR(i / 3 —> 77) essentially equal to that of the benchmark hf model. However, 
in general H® and H ° mix through the following mass matrix written in the {HfH® ) 
basis |211 


/ 8c^(Ai + A 3 ) 2\^6shCh^3 \ 2 

v 2\/6sfrCtfA 3 3s 2 h (\2 + A 3 ) J 


(19) 


Here A* are dimensionless quartic couplings in the Higgs potential and sh = sin 6h,ch = 
cos Oh- The assumption that the A; couplings are roughly the same order of magnitude 
together with the imposition of the bound sh < 0.4 results in very small mixing |3U. 
Moreover, H® would be the lightest Higgs boson in the HTM limit of small . 577 , as stressed 
in [32j. In this paper we will study the production process qq' —> H^H°', assuming that 
H± is a fcrmiophobic Higgs with BRs equivalent to the benchmark hf model. 


2.3 Fermiophobic Higgs boson branching ratios 

For the sake of illustration, we depict in Fig. |U the branching ratios of a fcrmiophobic 
Higgs boson hf into VV where V can be either a W, Z or 7. I 11 this figure we assumed 
that the hf couplings to fermions are absent and that hf ^ 77 is mediated solely by a 
W boson loop, 


h f 



7 


7 




7 


7 


giving rise to the following hf branching ratio into two photons, 


F {h f 77) 


a 2 g 2 

10247T 3 rriw 


Fi cos/3 | 2 


( 20 ) 


with Fi = F\(t), r = 4a function given in |2]. We remind the reader that the 
hfWW coupling normalized to the SM (p 0 WW coupling satisfies sin(/3 — a) —> cos/3 in 
the fermiophobic limit. 

This gives rise to benchmark BRs which are used in the ongoing searches to derive mass 
limits on rrif lf . In practice, hf —>77 can also be mediated by charged scalar loops: H ± in 
the 2HDM jTUj. XQ and Hf, Hf, Hf± in the HTM |33| . Although such contributions are 
suppressed relative to the W loops by a phase space factor, they can be important if the 
mixing angle suppression for the hfWW coupling (cos/3) is quite severe i.e. the scenario 
of interest in this paper. I 11 our numerical analysis we will assume the benchmark BRs 
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m hf [GeV] 


Figure 1: Branching ratios of the largest decay modes of a fcrmiophobic Higgs boson 
assuming exact fcrmiophobia at tree-level. The branching ratio into 77 equals the W*W* 
mode for rrif lf ps 90 GeV and drops to 20% for m/ l/ = 100 GeV. 

given in Fig. |T| One can see from the figure that the loop induced decay mode h/ —> 77 
is dominant for mh f ~ 95 GeV and drops below 0.1% for hf masses above 150 GeV. On 
the other hand, the decay channel hf —> W*W* dominates for rrih f g/95 GeV, being close 
to 100% until the threshold for hf decay into two real Z 1 s is reached. 

2.4 The decay H ± —> hfW* 

The experimental signature of the process qq' —> H ± hf depends on the decay modes of 
H ± . If H ± decays to two fermions then the signal would be of the type 77 + X, which 
is essentially the same as that assumed in the inclusive searches. However, crucial to our 
analysis is the fact that the decay H ± —> hfW* may have a very large BR [23’ in the 
2HDM (Model I). This is because the decay width to the fermions ( H ± —> /'/) scales 
as 1/tan 2 /3. Similar behaviour occurs in the HTM [25] for the decay Hf —> HfW* 
with the replacement l/tan 2 /3 —> tan 2 9h- Thus in the region of tan/3 > 10 (or small 
sin Q 11 ) the fermionic decays of H ± are depleted. This enables the decay H ± —> hfW* 
(Hf — > H\ y W*) to become the dominant channel even if the mass difference m#± — rrih f 
is much less than mw ■ In Fig. [2]we show the branching ratios of the charged Higgs boson 
into fermions and hfW* as a function of M H ± for several values of tan/3 and rrif lf . From 
the right panels we see that in the large tan/3 regime the fermionic decays are indeed 
suppressed. Moreover, we also see that for light fermiophobic Higgs bosons, where a W 
boson can be produced on its mass shell, the decay H ± —> W ± hf is essentially 100% 
for any tan /3. On the other hand, for heavier fcrmiophobic Higgs bosons, the fermionic 
decays can be the preferred decay channels mainly for small tan /3. 
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Figure 2: The charged Higgs boson branching ratios into fermions tv t (green/dotted), tb 
(red/dashed), cs (blue/dot-dashed), and cb (magenta/solid), and W*hf (black/solid) as a 
function of the charged Higgs boson mass for two different tan /3 and fermiophobic Higgs 
mass values. 

3 Phenomenology of hf at hadron colliders 
3.1 hf production via the VVhf coupling 

Current searches at the Tevatron assume that production of hf proceeds via the VVhf 
coupling (V = W, Z) that originates from the kinetic part of the Lagrangian. Run I 
searches utilized the process qq' —> Vhf giving a signature of a 77 and a vector boson 
ra.int The preliminary Run II search is for inclusive 77 2111] and is therefore sensitive 
to both qq' —» Vhf and the subdominant vector boson fusion qq' —► hfqq] see [TBJ. Note 
that gg —► hf via a fermion loop does not contribute to hf production. 

In the 2HDM (Model I) the strength of the VVhf coupling relative to the SM coupling 
VV(j)° is given by 

VVh f oc 1 . (21) 

y 1 + tan 2 f3 

Hence the production mechanism qq —> V* —► Vhf can be rendered completely ineffective 
for tan f3 > 10. In the HTM the fermiophobic Higgs boson has a coupling size relative to 
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VV(f>° given by 

WH? (22) 

In direct analogy to the large tan /3 case of the 2HDM (Model I), a small sh would suppress 
the coupling VVH® and consequently deplete the hfV production. Hence it is of concern 
to consider other production mechanisms which are unsuppressed in the above scenario. 


3.2 Associated hf production with H ± and the multi-photon sig¬ 
nature 


The production mechanism qq' —> H ± hf is complementary to that of qq' —> Vhf. This 
can be seen immediately from the explicit expressions for the couplings. In the 2HDM 
(Model I) one has 


while in the HTM 


VH^hf oc 


r± u 0' 


tan/3 

\Jl + tan 2 /3 

(23) 

2V2 

K y/3 CH ' 

(24) 


Hence the above couplings are unsuppressed in the region of the parameter space where 
the standard production mechanism qq' —> Vhf becomes ineffective. The larger coefficient 
for the VHfH® coupling is a consequence of the quantum number ( I,Y ) assignments in 
the HTM. 

To date complementary mechanisms have not been considered in the direct fermiopho- 
bic Higgs searches at the Tevatron. As emphasized in |22), |2H| a more complete search 
strategy for hf at hadron colliders must include such production processes in order to 
probe the scenario of fermiophobic Higgs bosons with a suppressed coupling hfVV. In 
the HTM one expects to be the lightest Higgs boson for small sin djj, which further 
motivates a search in the complementary channel qq' —> . 

The experimental signature arising from the complementary mechanism qq' —> H ± hf 
depends on the H ± decay channel. In a large fraction of the parameter space where the 
complementary mechanism qq' —> H±hf is important, the H ± decay is dominated by 
H ± —> hfW*. Consequently, this scenario would give rise to double hf production, with 
subsequent decay of hfhf 7777 , W 77 and VVVV. For light hf ( rrih f ~ 90 GeV), the 
signal 7777 would dominate, as discussed in [2IJ at LEP, in [22] for the Tevatron Run II 
and ^ at the LHC and a Linear Collider. More specifically, the multi-photon signature 
arises in the portion of the parameters space where m/ l/ 1 $ 90 GeV, m H ± 200 GeV, and 
tan/3 > 1 in the 2HDM Model I framework. In that region, BR(hf —> 77 ) ~ 1 and 
BR^H* hfW**) ~ 1 as well, leading to a 47 + leptons or jets signature. 

As explained in 123, processes other than qq' —> H ± hf could give rise to a 47 + X 
signal. One such mechanism is qq —> A°hf, where A 0 is the heavy neutral pseudoscalar 
decaying A 0 —> hfZ*. However, LEP already searched for e + e _ —> hfA° and set the bound 
rrih f + A > 160 GeV [H[. Thus any contribution from qq —> A°hf will be phase space 
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suppressed relative to that originating from qq' —> H ± hf. A similar argument applies to 
the production of a pair of charged Higgs bosons and its subsequent decay into hfV* pairs, 
i.e. qq —> Z*,^* —> H + H~ —> hfhfW + W~ which is phase space suppressed at Tevatron 
energies (2m#± > 180 GeV from direct searches). In the minimal supersymmetric 
model the total rates for H + H~ production are enhanced in the large tan /3 regime through 
the Yukawa couplings of Higgs bosons to bottom quarks [21], however in the 2HDM 
Model I and HTM, these Yukawa couplings are suppressed. The LHC would probably 
have much better prospects in these additional channels if all the above pair production 
mechanisms were combined with the H ± hf associated production in a fully inclusive multi¬ 
photon search. Since our analysis for the Tevatron we will focus on qq' —> H ± hf , which 
provides the best search potential for the very light hj region because the phase space 
constraint (m#± + rrih f > 100 GeV) is the least restrictive of the Higgs pair production 
mechanisms. 


4 Multi-photon signal analyses 

We now present our analysis for the inclusive production of multi-photon final states 
which may or may not be accompanied by extra leptons and/or jets, i.e. the reaction 

pp —> hfH ± —> hfhfW ± —> 7777 + X 

at the Tevatron Run II. We focus our attention on two inclusive final states; i) at least 
three photons (> 37 ) and ii) four photons ( 47 ). Only the “1-prong” tau lepton decays 
were considered. 

In our analysis we evaluated the signal and standard model backgrounds at the parton 
level. We calculated the full matrix elements using the helicity formalism with the help 
of Madevent 125]. We employed CTEQ6L1 parton distributions functions J2B] evaluated 
at the factorization scale Qf = where a/§ is the partonic center-of-mass energy. 
Although QCD corrections increase the tree-level cross section by a factor of around 
1.3 pa, we shall present results using the tree-level cross sections only. Moreover we 
included momenta smearing effects as given in [ISj, [ 2 E] and a detection efficiency of 85% 
per photon. 

We present in Fig. [21 the total signal cross section times the branching ratios of H ± —> 
hfW ± and hf 77 for the complementary process qq' —> H ± hf ; these results were 
obtained without cuts and detection efficiencies. In the left (right) panel we present the 4 7 
production cross section before cuts as a function of rrii lf for three different values of m#± 
and tan (3 = 3 (30). In the left panel, where tan/? = 3, the upper curve ( m H ± = 100 GeV) 
shows the strongest effect of the phase space suppression of the decay H ± —> W*hf for 
nih f ~ 60 GeV. This cross section reduction can be partially compensated by the increase 
in tan f3 as shown in the right panel. From the figure it is evident that this process will 
produce a large number of events before cuts over a large fraction of the parameter space. 

Potential SM backgrounds for the multi-photon signature of fcrmiophobic Higgs bosons 
are: i) the three and four photon production pp —> 777(7), ii) three photons and a W 
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Figure 3: Total production cross sections times branching ratios of hf —> 77 and H ± —► 
W ± hf for pp —> hfH ± —*■ hfhf + W ± 7777 + W ± before cuts at the Tevatron Run 11 
in femtobarns. The values of tan/? and m#± are as indicated in the figure. 

production pp —> 777 W, and iii) the associated production of two or three photons and a 
jet where the latter is misidentified as a photon pp —► 77(7 )j{—* 7 + -X"). We verified that 
after cuts and taking into account a P(j —> 7 ) = 4 x 10 -4 |38j photon misidentihcation 
probability the total SM background amounts to 3.8 events for an integrated luminosity 
of 2 fb -1 . Therefore the complementary process for the fermiophobic Higgs search has 
the great advantage of being extremely clean for a large portion of the 2HDM and HTM 
parameter space. In contrast, the ongoing search for inclusive 77 + X |2Tj. suffers from a 
sizeable background originating from QCD jets faking photons. For the exclusive channel 
(77 + V) the background is considerably smaller but still not negligible [TKj.jTTj. 

4.1 Searches at the Tevatron Run II 

The multi-photon topology is privileged concerning the level of background, which is small 
in the SM after mild cuts. Consequently, we imposed a minimum set of cuts on the final 
state particles, in order to guarantee their identification and isolation. Further studies 
could optimize the search strategy. We required the events to possess central photons 
with enough transverse energy to assure their proper identification 

> 15 GeV , [rf\ < 1.0 , (25) 

and isolated from the other particles in the final state (.X = charged lepton or jet) with 
a transverse energy in excess of 5 GeV 

> 0.4 , XR^x ^ 0.4 . (26) 

Notice that the high p? central photons is enough to guarantee the trigger of these events 
|38j . These cuts are very effective against the backgrounds from continuous 777 + X 
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Figure 4: Normalized transverse energy distributions (in GeV) of photons for two different 
charged Higgs boson masses and three different fermiophobic Higgs masses. The vertical 
solid lines indicate the Ej, cut in eq. (ESD- 


production which occur mainly through photon and gluon bremsstrahlung emission from 
initial and/or final state quarks, and gluon splitting to collinear quarks. We have checked 
that 4y + X topologies give a negligible contribution after imposing the cuts. 

In order to understand the effect of these cuts on the signal we studied some kinemat- 
ical distributions. We present in Fig. 0 the normalized transverse energy distribution of 
the final state photons for several values of Higgs masses and tan (3 = 30. As one can see, 
the Ej* spectrum peaks around < m,h f /2 and the spectrum at low E? decreases as the 
fermiophobic Higgs becomes heavier. We can also learn that the Ej, distribution becomes 
harder as the charged Higgs mass increases. Thus, the transverse energy cut in eq. (|251) 
attenuates more the light fermiophobic Higgs signal. 

Figure 0 contains the photon rapidity distribution for the same parameters used in 
Fig. Q] The rapidity distribution of the photons stemming from the fermiophobic Higgs 
decay peaks around zero. However, there is a sizeable contribution from high rapidity 
photons. For heavier charged Higgs bosons the rapidity distribution is more central. The 
hardest cut that we applied is the requirement that the absolute value of the photon 
rapidity be smaller than unity, and its effect is rather insensitive to the neutral Higgs 
mass. On the other hand, the separation cuts in eq. (HU) have little effect on the signal 
cross section as shown in Fig. |H1 with perhaps the exception of very small fermiophobic 
Higgs masses. Notice that we did not introduce any cut on the photon-photon invariant 
masses. Certainly if a signal is observed the photon pair invariant mass will display a clear 
peak at rrih f even after adding all the possible photon pair combinations and backgrounds; 
see Fig. 0 

We display in Figure |S] the region in the plane m#± 0 rrif lf where at least a 3cr signal 
can be observed, exhibiting three or more photons, in the framework of the 2HDM Model 
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Figure 5: Normalized rapidity distributions of photons for two different charged Higgs 
boson masses and three different fermiophobic Higgs masses. The vertical solid lines 
indicate the ?? 7 cut of eq. (1251) . 


1 for an integrated luminosity of 2 fb _1 . Statistical significance cr is defined by a = 
S/y/B, where S(B) is the number of signal (background) events after applying cuts and 
efficiency factors. A few comments are in order. First of all, the expected number of 
events diminishes for small hf masses since fewer events pass the Ej, cut in eq. ( 1251 ) as can 
be seen from Fig. 01 Secondly, the shape of the region presenting at least 3a (5a or 10a) 
significance in the large hf mass region is the result of a competition between the phase 
space suppression of the cross section as m#± increases for fixed m/ l/ and the growth of 
the H ± —> W ± hf branching ratio; see Fig. [2j Furthermore, for a fixed number of events, 
the optimum reach in rriH± takes place for rrih f — 30-40 GeV. This is a consequence of 
the combined effects of cuts and phase space suppression as we have already discussed. 

As one can see from the upper panel in Fig. |H1 even in the low tan /3 region the reach of 
the Tevatron RUN If is quite impressive in this scenario. If no events were observed above 
the backgrounds at RUN II, a large fraction of the m H ± ® m/ l/ plane would be excluded 
at the 3a level. The situation improves slightly for larger tan (3 as can be seen from the 
lower panel of Fig. |HJ Importantly, the expected number of events is rather large in the 
region rrih f < 70 GeV and m H ± < 150 GeV, such that it will be possible to reconstruct 
the hf mass from the photon-photon invariant mass distribution; see Fig. [71 

For comparison, we present in Fig. |51 the expected signal significance of events con¬ 
taining three or more photons after cuts for the HTM. In our numerical analysis we take 
ch = 1 as a benchmark value, and the signal significance for other values of ch can be 
obtained by simply rescaling the displayed numbers. From the bound Sh < 0.4 one ob¬ 
tains ch > 0.9. In the exact ch = 1 limit (i.e. triplet vev b = 0) the neutrinos would 
not receive a mass at tree-level (see eq. mi). Extremely small sh < 10 9 would require 
non-perturbative values of h t] to generate realistic neutrino masses. We are interested in 
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Figure 6: Normalized cone variable distributions between two photons (upper panels), a 
photon and a charged particle or jet (lower panels) for the same parameter used in Fig. 
0 The vertical solid lines indicate the A R cuts in eq. (I2H1) . 


the interval 0.9 < c# < 0.99 (corresponding to GeV scale triplet vev) in which H^ decays 
primarily to photons in the detector, and neutrino mass is generated with a very small 
hij ~ 10 - 10 . 

It is clear that a larger region of the rri) lf 0 m H ± parameter space can be probed in 
the HTM than in the 2HDM. In fact, at the 3cr level RUN II will be able to exclude 
Higgs masses up to m H ± < 240 GeV or m/ l/ 100 GeV. In order to understand the signal 
suppression if one requires an inclusive state containing four photons to pass our cuts, 
we present in Fig. El lower panel the expected number of events for the 2HDM, assuming 
tan/? = 30 and an integrated luminosity of 2 fb _1 . As expected, not only the reach in 
gets reduced to m H ± < 150 GeV at 95% C.L., but also the low and high hf mass 
regions become substantially depleted. 
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5 Conclusions 


Higgs bosons with very suppressed couplings to fermions (hf) can arise in various exten¬ 
sions of the Standard Model (SM) such as the Two Higgs Doublet Model (2HDM) Type 
I or Higgs Triplet Model (HTM). Their conventional production mechanism at hadron 
colliders qq' —> W ± hf can be severely suppressed by either large tan [3 or small triplet 
vacuum expectation value. In this scenario the complementary channel pp —> H ± hf is 
maximal and provides an alternative production mechanism. We studied the reaction 
qq' —> H ± hf followed by the potentially important decay H ± —* hfW*. We performed 
a Monte Carlo simulation of the detection prospects for a light hf where the branching 
ratio into photon pairs is dominant, which gives rise to multi-photon signatures with very 
low SM background. We showed that if a signal containing at least three photons is not 
seen at the Tevatron RUN II then a large portion of the m H ± versus rrif lf plane can be 
excluded both in the small and large tan /3 regimes of the 2HDM. Conversely, if a signal 
were observed then > 50 events are expected for a light H ± and hf, which would allow 
further detailed phenomenological studies. 


Acknowledgements 

We would like to thank Oleksiy Atramentov for discussions concerning the SM back¬ 
grounds. This research was supported in part by Fundagao de Amparo a Pesquisa do 
Estado de Sao Paulo (FAPESP), by Conselho Nacional de Desenvolvimento Cientffico e 
Tecnologico (CNPq), and by Conicyt grant No. 1040384. 


References 

[1] T. J. Weiler, Proceedings of the 8th Vanderbilt Int. Conf. on High Energy Physics, 
Nashville, TN, Oct 8-10, 1987: ; Edited by J. Brau and R. Panvini (World Scientific, 
Singapore, 1988), p219 

[2] H. E. Haber, G. L. Kane and T. Sterling, Nucl. Phys. B 161, 493 (1979) 

[3] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, “The Higgs Hunter’s Guide,” 
(Reading, MA: Addison-Wesley, 1989) 

[4] H. Georgi and M. Machacek, Nucl. Phys. B 262, 463 (1985); M. S. Chanowitz and 
M. Golden, Phys. Lett. B 165, 105 (1985) 

[5] A. Stange, W. J. Marciano and S. Willcnbrock, Phys. Rev. D 49, 1354 (1994) 

[6] M. A. Diaz and T. J. Weiler, arXivdiep-ph/9401259 

[7] V. D. Barger, N. G. Deshpande, J. L. Hewett and T. G. Rizzo, arXivdiep-ph/9211234. 
In Argonne 1993, Physics at current accelerators and supercolliders* 437-442 


16 



[8] H. Pois, T. J. Weiler and T. C. Yuan, Phys. Rev. D 47, 3886 (1993) 

[9] A. G. Akeroyd, Phys. Lett. B 368, 89 (1996) 

[10] A. Barroso, L. Brucher and R. Santos, Phys. Rev. D 60, 035005 (1999) 

[11] L. Brucher and R. Santos, Eur. Phys. J. C 12, 87 (2000) 

[12] G. Abbiendi et al. [OPAL Collaboration], Phys. Lett. B 544, 44 (2002) 

[13] P. Abreu et al. [DELPHI Collaboration], Phys. Lett. B 507, 89 (2001); Eur. Phys. 
J. C 35, 313 (2004) 

[14] A. Heister et al. [ALEPH Collaboration], Phys. Lett. B 544, 16 (2002) 

[15] P. Achard et al. [L3 Collaboration], Phys. Lett. B 534, 28 (2002); Phys. Lett. B 568, 
191 (2003). 

[16] B. Abbott et al. [D0 Collaboration], Phys. Rev. Lett. 82, 2244 (1999) 

[17] T. Affolder et al. [CDF Collaboration], Phys. Rev. D 64, 092002 (2001) 

[18] S. Mrenna and J. Wells, Phys. Rev. D 63, 015006 (2001) 

[19] V. M. Abazov et al. [DO Collaboration], arXiv:hep-ex/0508054 

[20] G. Landsberg and K. T. Matchev, Phys. Rev. D 62, 035004 (2000) 

[21] A. Melnitchouk [DO Collaboration], Int. J. Mod. Phys. A 20, 3305 (2005) 

[22] A. G. Akeroyd and M. A. Diaz, Phys. Rev. D 67, 095007 (2003) 

[23] A. G. Akeroyd, M. A. Diaz and F. J. Pacheco, Phys. Rev. D 70, 075002 (2004) 

[24] A. G. Akeroyd, Nucl. Phys. B 544, 557 (1999) 

[25] A. G. Akeroyd, Phys. Lett. B 442, 335 (1998) 

[26] J. F. Gunion and H. E. Haber, arXiv:hep-ph/0506227 

[27] J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980). 

[28] F. Cuypers and S. Davidson, Eur. Phys. J. C 2, 503 (1998); Y. Kuno and Y. Okada, 
Rev. Mod. Phys. 73, 151 (2001) 

[29] J. F. Gunion, R. Vega and J. Wudka, Phys. Rev. D 42, 1673 (1990) 

[30] A. Kundu and B. Mukhopadhyaya, Int. J. Mod. Phys. All, 5221 (1996); H. E. Haber 
and H. E. Logan, Phys. Rev. D 62, 015011 (2000) 

[31] A. G. Akeroyd, Phys. Lett. B 353, 519 (1995). 


17 



[32] P. Bamert and Z. Kunszt, Phys. Lett. B 306, 335 (1993) 

[33] J. F. Gunion, R. Vega and J. Wudka, Phys. Rev. D 43, 2322 (1991). 

[34] A. Alves and T. Plehn, Phys. Rev. D 71, 115014 (2005). 

[35] F. Maltoni and T. Stelzer, JHEP 0302, 027 (2003); T. Stelzer and W.F. Long, Phys. 
Commun. 81, 357 (1994). 

[36] J. Pnmplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky, W.K. Tung, JHEP 0207, 

012 ( 2002 ). 

[37] S. Dawson, S. Dittmaier and M. Spira, Phys. Rev. D 58, 115012 (1998) 

[38] A. Gajjar [CDF Collaboration], arXiv:hep-ex/0505046; V. M. Abazov et al. [DO 
Collaboration], Phys. Rev. D 71, 091108 (2005) A. Melnitchouk, Ph.D. Thesis, 
FERMILAB-THESIS-2003-23 (2003). 


18 






m yy [GeV] myy [GeV] 


Figure 7: The photon-photon invariant mass spectrum for rn h± = 150 GeV, rri} lf = 30, 
50, 70 GeV, and m H ± = 140 GeV, m/ l/ = 90 GeV at the upper left, upper right, bottom 
left, and bottom right panels, respectively. Also shown are the sum of all backgrounds 
consisting of 3y + X final states. In these plot we entered the invariant mass of all photon 
pair possible combinations. In all cases we set tan [3 = 30. 
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Figure 8: Expected signal statistical significance presenting three of more photons in the 
rriH± ® mh f plane assuming an integrated luminosity of 2 fb -1 at the Tevatron RUN II. 
We assumed the 2HDM Model I and took tan/? = 3 (30) in the upper (lower) panel. 
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Figure 9: In the upper panel we display the expected signal statistical significance con¬ 
taining three of more photons in the m h± ®mh f plane in the HTM framework, assuming 
an integrated luminosity of 2 fb -1 at the Tevatron RUN II. In the lower panel we present 
the expected number of events presenting four photons in the m#± ® plane for the 
2HDM Model I and assuming tan/3 = 30. 
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